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We study peculiarities of proximity effect in clean superconductor - ferromagnet structures caused 
by either spatial or momentum dependence of the exchange field. Even a small modulation of the 
exchange field along the quasiparticle trajectories is shown to provide a long range contribution 
to the supercurrent due to the specific interference of particle- and hole- like wave functions. The 
momentum dependence of the exchange field caused by the spin - orbit interaction results in the 
long - range superconducting correlations even in the absence of ferromagnetic domain structure 
and can explain the recent experiments on ferromagnetic nanowires. 



The exchange field h in ferromagnetic (F) metals is 
well known to destroy Cooper pairs resulting, thus, in a 
strong decay of superconducting (S) correlations in the 
F material and suppression of Josephson current in SFS 
junctions (see Refs. [J, for review). Considering the 
quantum mechanics of quasiparticle excitations this de- 
structive effect of the exchange field can be viewed as a 
consequence of a phase difference 7 ~ L/£h — 2Lh/hVF 
gained between the electron- and hole- like parts of the 
total wave function at the path of the length L. Both 
in the clean and dirty limits the measurable quantities 
should be calculated as superpositions of fast oscillating 
contributions e* 7 from different trajectories and, thus, 
rapidly vanish with the increasing distance from the SF 
boundary. 

This textbook physical picture appears to be in sharp 
contrast with a number of recent experiments [3 8] which 
point to an anomalously large length of decay of super- 
conducting correlations inside the F metal. As we can 
judge from the observation Q of a noticeable supercur- 
rent through a Co nanowire, this decay length can be of 
the order of half a micrometer which well exceeds typi- 
cal coherence lengths in ferromagnets both in the clean 
and dirty limits. In the dirty limit such strong proxim- 
ity effect can hardly be explained even taking account 
of long-range triplet correlations Q induced by the ex- 
change field inhomogeneity. 

Naturally, the inhomogeneity of the field h caused 
by the ferromagnetic domain structure can improve the 
conditions of Cooper pair survival in the clean limit as 
well. To suppress the destructive trajectory interference 
mentioned above the domain structure should cancel the 
phase gain 7 for a certain group of quasiparticle trajec- 
tories. A simple example of such phase gain compen- 
sation can be realized in a clean junction consisting of 
two F layers with opposite orientations of magnetic mo- 
ment 0, [l(|. On the other hand in the diffusive limit 
this compensation effect vanishes Note, that the ex- 
change field inhomogeneity along the quasiclassical tra- 



jectory experiencing multiple reflections from the ferro- 
magnet surface can appear even in the absence of the 
spatial domain structure. Indeed, the exchange field act- 
ing on band electrons in a solid with a finite spin - orbit 
interaction should obviously depend on the quasiparticle 
momentum [12j: h = h(k). The normal quasiparticle 
reflection is accompanied, of course, by the change in the 
momentum direction, and, thus, by the change in the 
exchange field. The momentum dependent h field can 
strongly affect the phase gain 7 along the trajectories 
even in the F sample prepared in a single domain state 
(as it has been done in experiments with Co nanowires 
I). 

The goal of this paper is to show that in the clean 
limit there exists a possibility to cancel the particle - 
hole phase difference for a large group of quasiclassical 
trajectories due to either spatial or momentum depen- 
dence of the exchange field. Such set of trajectories pro- 
vides a long-range contribution to the Josephson cur- 
rent through a ferromagnetic system which decays at the 
length scale characteristic for a nonmagnetic metal. We 
consider two generic examples which illustrate the above 
scenario of a long-range proximity effect: (i) Josephson 
transport through a pair of ferromagnetic layers with a 
stepwise exchange field distribution; (ii) Josephson trans- 
port through a nanowire with a specular electron re- 
flection at the surface and exchange field varying with 
the changing quasiparticle momentum. See Supplemen- 
tal Material at [URL will be inserted by publisher] for 
details of calculations. 

Josephson transport through a ferromagnetic bilayer. Let 
us start from the simplest model illustrating the origin 
of the quasiparticle interference suppression: Josephson 
junction containing two ferromagnetic layers of the thick- 
nesses di and d 2 , respectively (see Fig.l). Here we con- 
sider the limit of short junction di + ffe <C £ s , where £ s 
is the superconducting coherence length. The exchange 
fields hi and h 2 in the layers are rotated at the angle 
a. Following the quasiclassical procedure considered in 
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Figure 1: (Color online) Josephson junction containing two 
ferromagnetic layers. Linear quasiparticle trajectory is shown 
by the red dashed line. 



Ref. 13) we find the current - phase relation: 

((n, np ) cos Try) 
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where n is the unit vector normal to the junction plane, 
rip is the unit vector along the trajectory, and a n are 
the coefficients of the Fourier expansion for the current 
- phase relation IsNsi'fi) for zero exchange field, i.e., 
for superconductor - normal metal junction of the same 
geometry. The angular brackets denote the averaging 
over different quasiclassical trajectories. The first two 
coefficients in this expansion take the form: 
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Ao is the 
gap, N = 

= (W2tt) 2 ) 
. ds is taken 



where fi m = 27r 2 T 2 (2m + 1) 2 /A^ + 1 
temperature dependent superconducting 
Sq 1 f ds J dnp(np, n), s^ 1 = kp/2n (sq 1 
for 2D (3D) junctions, and the integral J . 
over the junction cross-section. The factor N is deter- 
mined by the number of transverse modes in the junction: 
N ~ S/so, where S is the junction cross-section area. 

The phase 7 can be found from the singlet part of the 
anomalous quasiclassical Green function: 

f s (s = s R ) = cos 7 

taken at the right superconduc ting electrode. Here we 
use a standard parametrization LJ] / = f s + f t a, where 
<7 is a Pauli matrix vector in the spin space. The functions 
f s , ft satisfy the linearized Eilenberger equations written 
for zero Matsubara frequencies 

- ihV F d s f s + 2hf t = , -ihV F d s { t + 2f s h = , (3) 

and the conditions f s (s = sl) — 1, tt(s = sl) = at 
the left superconducting electrode. Solving the above 
equations for the particular bilayer geometry we find: 
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where cos0 — (n,nj?). This expression allows us to write 
the first harmonic in the current - phase relation in the 
form: 
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where I c i(d/Ch) is the critical current of the first har- 
monic in a SFS junction with a homogeneous exchange 
field h. The interference effects discussed in introduction 
result in the power decay of the critical current I cl vs 
the F layer thickness d: I c \ oc d" 1 / 2 for a 2D junction 
15] and I c i cx d^ 1 for a 3D junction [lfjj]. Taking sym- 
metric case d\ = di we immediately get a long-range 
contribution to the Josephson current 



Sid = sin 2 — Id (0) sin ip 



(6) 



which does not decay with the increasing distance be- 
tween the S electrodes. It is important to note that this 
contribution does not vanish for an arbitrary nonzero an- 
gle between the magnetic moments in the F layers. 

Long-range behavior can be observed for a second har- 
monic in the current - phase relation as well. Indeed, 
calculating the average ((n, n^) cos 27) we find a non- 
vanishing long-range supercurrent contribution even for 
di ^ d 2 : 
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Note, that the emergence of long-range proximity effect 
for high harmonics in Josephson relation is in a good 
agreement with recent theoretical findings in Refs.ll7l.ll8. 
Josephson current through a ferromagnetic wire. We now 
proceed with the consideration of a more complicated 
example of the interference phase suppression in a fer- 
romagnetic wire where the quasiclassical trajectories of 
electrons and holes experience multiple specular reflec- 
tions from the wire surface (see Fig. 2a). The particular 
geometry shown in Fig. 2a can be considered as a rough 
model for experiments on Co nanowires @- For simplic- 
ity we restrict ourselves to the case of a 2D junction. 

Taking account of the spin-orbit interaction inside the 
ferromagnet we obtaine the exchange part of the effec- 
tive Hamiltonian for the band electrons depending on 
the quasi-momentum (k) orientation: 

Hex = J3 /^ (k^oiOj = h(k)a , 

ij 

where ho is a pseudo vector determined by the ferro- 
magnetic moment. Assuming the absence of the system 
anisotropy described by a polar vector we find the sim- 
plest form of the resulting exchange field: h ~ ho + 
Psokp 2 (ho, k)k, where (3 so is a constant determined by 
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orbit interaction, and fcp- is the Fermi momen- 



3 



D 



y 



a) 



s L 




s R 

l D 

3 >X 






C 5 

L 





b) 



Figure 2: (Color online) Josephson transport through a 
nanowire in the overlap (a) and edge (b) geometries. The 
quasiparticle trajectories are shown by the red dashed lines. 



The exchange field along the quasiparticle trajectory 
experiencing the reflection at the wire surface should 
change its direction. Thus, we obtain the problem de- 
scribed by Eqs. © with a periodic exchange field along 
the trajectory characterized by a given angle 9 and a cer- 
tain starting point at the superconductor surface. The 
same equations for each trajectory can be of course de- 
rived for a periodic domain structure. Let us consider 
first the problem of calculating the band spectrum e(k) 
in the field h varying with the period 2Dj sin#: 



-ihV F d s f s 
-ihV F d s i t 



-2M t = e{k)f s , 
2/ s h - e(k)f t . 

The solution can be written in the Bloch form: 
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where f sk (s + 2D/ sm.6) = f s k(s) and f tk {s + 2D/ sm9) = 
f*fe(s). One can see that provided this solution corre- 
sponds to the energy branch e a (k) there exist another so- 
lution (/*, — f£) corresponding to the energy — e a (k). On 
the other hand the latter solution corresponds also to the 
energy es-(-k) and, thus, we obtain the following sym- 
metry property of the band spectrum: eg-(— fc) — —e a (k), 
where the indices a and a denote different branch num- 
bers. The full set of energy branches can be split in 
such pairs provided the number of branches is even. For 
an odd number of branches there is always one branch 
which does not have a partner. For this branch we ob- 
tain e a {— k) = —e a (k) and, thus, this spectrum branch 
crosses the zero energy level at k — 0: e CT (0) = 0. The 
corresponding phase gain 7 appears to vanish for trajec- 
tories containing an integer number of periods shown in 
Fig. 2a and, therefore, the solution with k = and e = 
provides a long-range contribution to the supercurrent. 

For the sake of definiteness we choose the field ho to be 
directed along the wire axis x and obtain the exchange 



field in the form: h = Ttoh x +yoh y (s), where h x (9) ~ ho 
is constant along the trajectory and h y (s) is a periodic 
function with zero average. In the interval —D/sm9 < 
s < D/ sin# the h y field component is defined by the 
expression h y = /3 so ho sin 9 cos signs. Introducing the 
Fourier expansions 

we rewrite the Eqs. © and © in the form: 



(hV F (k + q)- e)F s (k + q) + 2h x F x (k + q) 

+2^H g _ i} F y (k + q) = 0, (10) 

(hV F (k + q)- e)F x (k + q) + 2h x F s (k + q) = , (11) 
(hV F (k + q) -e)F y {k + q) 

+2^H q _ s F s (k + q)=0. (12) 

Here q, q — q m — 7r(2m+ 1) s'm9/ D, m is an integer, and 
h — j3 so ho sin 9 cos 9. 

To get the solution for a small periodic field h y we 
use a perturbative approach similar to the nearly free 
electron approximation in the band theory of solids and 
restrict the number of interacting Fourier harmonics in 
the expansions. For this purpose it is instructive to con- 
sider the limit of zero periodic potential h y and separate 
three types of solutions: (i) the solution (F s ,F x ,F y ) = 
(0, 0, l)5 q - p corresponding to the energy eo = hV F (k + p) 
(ii) the solutions (F s ,F x ,F y ) = (1, ±1, 0)6 q - p± corre- 
sponding to the energies e± = hV F (k +p±) ± 2h x . Here 
p and p± are arbitrary reciprocal lattice vectors. The 
above modes should strongly interact provided the res- 
onant condition eo = e + = e - is fulfilled. Such reso- 
nance is possible for the case when the value 2h x /hV F 
equals to a certain reciprocal lattice vector q m . Close 
to such Bragg - type resonance we see that the dom- 
inant harmonics correspond to the following choice of 
reciprocal lattice vectors: p = 0, p± = T'Zm- Writing 
the solution as a superposition of these three harmon- 
ics we find renormalized spectral branches eo = HV F k, 
e± = HV F k ± yJ{hV F q m — 2h x ) 2 + 8|7? 9m | 2 and corre- 
sponding eigenfunctions. Applying now the boundary 
conditions at s = for the superposition of the above 
eigenfunctions we find the amplitude of the singlet com- 
ponent corresponding to the energy branch eo and k = 0: 
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At the surface of a right superconducting electrode we 
should take the coordinate s to be equal to the integer 
number of periods. We also need to sum up the above 
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resonant expressions over all Fourier harmonics of the 
periodic potential: 



fs(s = S R )=y2 ~7TTr 
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The precision of such resonant - type expression has 
been also confirmed by the numerical solution of the 
Eqs. © and ^ carried out using the transfer matrix 
method. Note, that we omit here the contribution from 
the solutions corresponding to the branches e±: these 
functions correspond to a nonzero quasimomentum and, 
thus, should gain a finite phase factor along the trajec- 
tory length. During averaging over different trajectories 
this phase factor causes the suppression of the resulting 
supercurrent contribution with the increasing wire length 
L. 

The starting point of the trajectory varies in the in- 
terval Ax = 2D f tan 6* and, as a consequence, the long - 
range first harmonic in current - phase relation takes the 
form: 

7T/2 

I\ = a\ sin if J d6 cos 6f s (sr) . 
a 

Assuming the resonances to be rather narrow we approx- 
imate them by the delta - functions and obtain: 

. ^ y/2nhV F h(e m ) . a 
h = ax sin if j^jj sin m . 

m x 

where sin^m = 2h x D / VSVp^m +1). In the limit D 
HVp/2h x one can replace the sum over m by the integral: 
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Certainly, the above long-range effect in the first har- 
monic is rather sensitive to the system geometry: taking, 
e.g., the system sketched in Fig. 2b we will not obtain the 
full cancellation of the phase 7 because the trajectories 
in this case do not contain integer number of exchange 
field modulation periods. However, similarly to the case 
of bilayer the long-range effect is still possible for higher 
harmonics. We apply the above perturbative procedure 
for the calculation of the full f s function for the geometry 
shown in Fig. 2b. The second harmonic in the current- 
phase relation reads 

tt/2 

7 2 = a2 sin2^ J decose(2{f 2 s { SR )) yo ~l) , (13) 



where (. . .) yo = (1/-D) ■ ■ ■ dyo denotes averaging over 
the starting point of the trajectory y (see Fig. 2b). 



Keeping only the terms linear in the small | H qm | ampli- 
tude we get the following expression for the long-range 
part of the second harmonic I2' 



I2 = a.2 sin 2ip 



\ - V2nhV F h(0 m ) , 2 

sin Tt 



hlD 



a2-^-Pso sin 2^. 



We emphasize that the second harmonic of Josephson 
current in both above examples is negative because of 
the condition a 2 < 0. 

Note that the absence of the decay of the sing le-channel 
critical current was pointed out in Ref. 19 as a pos- 
sible source of the long-ranged proximity effect in Co 
nanowires. However the averaging of the phase gain for 
different modes strongly decreases the critical current. 
In contrast the results presented in this Letter demon- 
strate that in the ballistic regime the spin-orbit inter- 
action generates the non-collinear exchange field which 
produces the long - range Josephson current. This con- 
clusion is always true for the second harmonic in the cur- 
rent - phase relation and for some geometries it may be 
also valid for the first harmonic. Therefore our findings 
provide a natural explanation of the recent experiments 
with Co nanowire @. To discriminate between two pro- 
posed mechanisms of the long ranged effect, the studies 
of higher harmonics in Josephson current-phase relations 
could be of major importance. Also it should be interest- 
ing to verify on experiment the predicted simple angular 
dependence (J6j> of the critical current in S/F/S junctions 
with composite interlayer. 
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